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Abstract 

Complex processes in perforated domains occur in many real-world applications. These problems 
are typically characterized by physical processes in domains with multiple scales (see Figure for the 
illustration of a perforated domain). Moreover, these problems are intrinsically multiscale and their 
discretizations can yield very large linear or nonlinear systems. In this paper, we investigate multiscale 
approaches that attempt to solve such problems on a coarse grid by constructing multiscale basis functions 
in each coarse grid, where the coarse grid can contain many perforations. In particular, we are interested 
in cases when there is no scale separation and the perforations can have different sizes. In this regard, we 
mention some earlier pioneering works [n i ir ^ ITT]. where the authors develop multiscale hnite element 
methods. In our paper, we follow Generalized Multiscale Finite Element Method (GMsFEM) and develop 
a multiscale procedure where we identify multiscale basis functions in each coarse block using snapshot 
space and local spectral problems. We show that with a few basis functions in each coarse block, one 
can accurately approximate the solution, where each coarse block can contain many small inclusions. 

We apply our general concept to (1) Laplace equation in perforated domain; (2) elasticity equation in 
perforated domain; and (3) Stokes equations in perforated domain. Numerical results are presented 
for these problems using two types of heterogeneous perforated domains. The analysis of the proposed 
methods will be presented elsewhere. 


1 Introduction 

Among multiscale problems, the problems in perforated domains are of great interest for many applications. 
The main characteristics of these problems is that the underlying processes occur in multiscale domains 
where the geometry of the domain has multiple scales, e.g., domain outside inclusions (see Figure [^. There 
are many important applications for processes in perforated domains. For example, fluid flow in porous 
media, diffusion in perforated domains, mechanical processes in hollow materials, and so on. Typically, it 
is a combination of a physical process and a heterogeneous media that gives rise to problems in perforated 
domains, e.g., fluid flow in porous media is a problem in perforated domain, while heat conduction in porous 
media may not be a problem in perforated domain. 
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The problems in perforated domains are of multiscale nature. The solution techniques for these problems 
require high resolution. In particular, the discretization needs to honor the irregular boundaries of perfora¬ 
tions. This gives rise to a hue-scale problems with many degrees of freedom which can be very expensive 
to solve. In this paper, our goal is to develop coarse-grid approaches where the coarse grids do not have to 
align with perforations and the perforated domains do not need to have a scale separation. Our objective is 
to construct a low dimensional approximate model for solving hne-scale problems. 

When the perforated domains have some scale separation, there are numerous works which include the 
works on the homogenization and asymptotic expansion in periodic perforated domains [T6l[3l|4]. Recently, a 
novel method of meso-scale asymptotic approximation for Laplace operator is introduced ([El) for the cases 
with a large number of perforations. The multiscale analysis of the perforated domain has been conducted 
using heterogeneous multiscale finite element method (HMMFEM) and Mulitiscale Finite Element Method 
(MsFEM) [El HH flSj - The authors in [m ng uni have extended multiscale finite element methods to 
arbitrary perforated domains and presented novel numerical approaches using Crouzeix-Raviart coupling of 
multiscale finite element basis. Their approach allows perforations to intersect the boundaries of coarse-grid 
blocks. In their approach, the authors use oversampling technique and weak coupling via Crouzeix-Raviart to 
avoid directly imposing boundary conditions for the basis functions. The latter can be difficult for multiscale 
basis function construction in perforated domains. In our paper, we also follow this general concept and 
avoid constructing boundary conditions for multiscale basis functions and construct them via local spectral 
decomposition. Our approach follows a Generalized Mulitiscale Finite Element Method (GMsFEM) [TT] in 
constructing multiscale basis functions which we discuss next. 

GMsFEM is a general multiscale procedure where multiscale basis functions in each coarse grid are 
constructed and is designed for many applications lainiE]. In the GMsFEM framework, one divides the 
computations into two stages, i.e., the offline stage and the online stage. In the offline stage, a reduced 
dimensional space is constructed, and it is then used in the online stage to construct multiscale basis functions. 
These multiscale basis functions can be re-used for any input parameter to solve the problem on a coarse 
grid. The main idea behind the construction of the offline and online spaces is to design appropriate snapshot 
spaces and determine an appropriate local spectral problem to select important modes in the snapshot space. 
In m, several general strategies for designing the local spectral procedures have been proposed. To apply 
GMsFEM for problems in perforated domains, one needs a novel way to construct the snapshot space and 
local spectral problem. In earlier works [Dllillig, heterogeneous problems have been studied. The main 
challenge is to take into account the domain heterogeneities when constructing multiscale basis functions. 
Gompared to problems in heterogeneous media without perforated regions, we impose boundary conditions 
on perforations. We generate local snapshot space using harmonic extensions of appropriately chosen a large 
set of boundary conditions. Local spectral problems are used to select dominant modes. We also discuss the 
use of randomized snapshots ([2]) that allow calculating a few snapshot vectors. The latter is important in 
reducing the computational cost. 

The above procedure is formulated for a general problem in perforated media. In this paper, we apply it 
to three different problems that arise in many applications. These include: (1) Laplace equation in perforated 
domain; (2) elasticity equation in perforated domain; and (3) Stokes equations in perforated domain. We 
show that our general concept can be applied to all these problems. We consider several numerical examples 
using two sets of perforated domains. One domain contains several large inclusions (with the sizes comparable 
to the coarse-grid sizes) and the other domain contains small inclusions (with the sizes much smaller than the 
size of the domain). We solve the above three PDFs in these two sets of domains and compare the multiscale 
solutions with the fine-grid solutions. Our preliminary numerical results demonstrate a fast convergence of 
the proposed GMsFEM approaches. In particular, we show that one can achieve a good accuracy with a 
very few degrees of freedom, especially for Laplace and elasticity problems. For Stokes equations, there is a 
room for improvement by enriching the pressure space. We discuss it in our numerical results. 

The rest of the paper is organized as follows. In Section we present preliminaries on the problems in 
perforated domains and the GMsFEM framework. The construction of the coarse spaces for the GMsFEM 
is discussed in Section [^ In Section numerical results for several representative examples are presented. 
Finally, we conclude our paper with some remarks in Section 
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2 Preliminaries 


In this section, we first present the underlying problem and the corresponding fine-scale discretization. Then, 
we discuss the multiscale strategy of solving this problem. We consider 


C^{u) = f in 

(1) 

3ij 

= 0 or -— = 0, on fl dB^ , 

on 

(2) 

u = g, on dVt D dVf . 

(3) 


Here, G (d = 2,3) is a bounded domain covered by inactive cells (for Stokes flow and Darcy flow) or 
active cells (for elasticity problem) . We call active cells where the underlying problem is solved, while 
inactive cells being the rest of the region. Suppose the distance between inactive cells (or active cells) is 
at most 5 and we use e scripts to denote the perforated domains. Denote the remaining part as i.e. 

= Q\B^. See Figure for an illustration of the perforated domain. denotes a linear differential 
operator, e.g. C^{u) = —Au for Darcy flow, n is the unit outward normal to the boundary and / G F{Q^) 
and g G G{Cl^) denote functions with a suitable regularity. Note that 5 is used to denote heterogeneous 
domains, even though we do not assume periodicity or scale separation. 



Figure 1: Illustration of a perforated domain. 

To simplify the notation, we denote by the appropriate solution space, and 

W(O^) = {^ G V{n^),v = 0 on dn^}. 

The variational formulation of Problem § is to find u G V{Q^) such that 

= if, for all v e 

(•, denotes a specific inner product over for either scalar functions or vector functions. In the following, 
we give some specific examples for the above abstract notations. 

• For the Laplace operator with homogeneous Dirichlet boundary conditions on 50^, we have 

C^{u) = —Au, (4) 

and H(D^) = {C^{u),v)^e = {Vu,Vv)^e- 
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• For the elasticity operator with homogeneous Dirichlet boundary condition on 9(1^, we let u G 
be the displacement field. The strain tensor e{u) G is defined by 

e{u) = + Vu^). 

In this paper, we assume the medium is isotropic. Thus, the stress tensor (j{u) G relates 

to the strain tensor s{u) in the following way 

a{u) = 2fie + AV • ul^ 

where A > 0 and /i > 0 are the Lame coefficients. We have 


jC^{u) = -V-(7, (5) 

where V{Q^) = and {C^{u)^v)^e = {(j{u)^cF{y))^s- 

• For Stokes equations, we have 


C^u ,p) 


fVp — jiiAu\ 

V v-« T 


(6) 


where /i is the viscosity, p is the fluid pressure, u represents the velocity, ^((7^) = x Lim, 

and 


{C^u ,p),{v ,q))^. 


/(Vu,Vv)n= 


We recall that Lq{Q,^) contains functions in with zero average in $7® 


We now turn our attention to a numerical approximation of the variational problem above. Due to the nature 
of this multiscale problem we consider the framework of multiscale finite element method m- The idea is 
computing the solution on a coarse grid instead of calculating on the fine grid directly. In the following, we 
introduce the necessary concepts and notations. 

Let be a coarse-grid partition of the domain and 7^ be a conforming fine triangulation of We 
assume that 77 is a refinement of 77/, where h and H represent the fine and coarse mesh sizes, respectively. 
Typically, we assume that 0 < h <C iL < 1, and that the fine-scale mesh 77 is sufficiently fine to fully resolve 
the small-scale information of the domain while i7 is a coarse mesh containing many fine-scale features. On 
the triangulation 77, we introduce the following finite element spaces 


Vh := {v e V{n^)\v\K e {P\K)y for all K € %}, 


where, k = 1, 2, denotes the polynomial approximation space, and 1 = 1,2 indicates either a scalar or a 
vector. 

In this paper. Generalized Multiscale Finite Element Method (GMsFEM) framework is applied. In the 
GMsEEM methodology, one divides the computations into offline and online computations. The offline 
computations are based on a preliminary dimension reduction of the fine-grid finite element spaces (that 
may include dealing with additionally important physical parameters, uncertainties and nonlinearities), and 
then the online procedure (if needed) is applied to construct a reduced order model. We start by constructing 
offline spaces. 

We construct the coarse function space 


Vos ■■= spanj^J^i, 

where N is the number of coarse basis functions. Each (pi is supported in some coarse neighborhood wi (see 
Eigures|^and[^for an illustration of coarse neighborhoods). 
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The overall performance of the resulting GMsFEM depends on the approximation properties of the 
resulting offline and online coarse spaces directly. Usually, a spectral problem is involved for a good approx¬ 
imation of the local solution space. In this paper, we focus on the construction of the offline spaces since the 
differential operators is parameter independent and linear. 

The GMsFEM seeks an approximation uq G which satisfies the coarse-scale offline formulation, 

{C^{uo),v)^e = {f,v)ns for all v e V°^. (7) 

Recall that the definitions of the bilinear forms {jC^(uq), v )are defined above, and (/, is the L 2 inner 
product. 

We can interpret the method in the following way using matrix representations. Recall that the coarse 
basis functions are defined on the fine grid, and can be represented by the fine-grid basis functions. 

Specifically, we introduce the following matrices: 

Rq = [01, • • • , ^ATe], 

where we identify the basis 0^ with their coefficient vectors on the fine-grid basis. Then, the matrix analogue 
of the system 0 can be equivalently written as 

= RoF, ( 8 ) 

where F is the fine-grid discretization of /. Further, once we solve the coarse system ([^, we can recover 
the fine scale solution by RqUq. In other words, Rq can be regarded as the transformation (also known as 
interpolation or downscaling) matrix from the space Uoff to the space 14 . 

The accuracy of the GMsFEM relies on the coarse basis functions {0i}. We shall present the construction 
of suitable basis functions for the differential operators in Section 

3 Local basis functions 

In this section we describe the offline-online computational procedure, and elaborate on some applicable 
choices for the associated bilinear forms to be used in the coarse space construction. Below, we offer a 
general outline for the procedure. 

1. Offline computations: 

— 1.0. Goarse grid generation. 

— 1.1. Gonstruction of a snapshot space that will be used to compute an offline space. 

— 1.2. Gonstruction of a small dimensional offline space by performing dimension reduction in the 
space of local snapshots. 

2. Online computations: 

— 2.1. For each input parameter, compute multiscale basis functions (it is needed for parameter- 
dependent or nonlinear problems). 

— 2.2. Solve a coarse-grid problem for a forcing term and boundary condition. 

— 2.3. Iterative solvers, if needed. 

In the offline computation, we first construct a snapshot space Unap or Ugnap^ depending on the choice of 
domain to generate the snapshot space, where cjd is an oversampled region that contains a coarse neighbor¬ 
hood uoi. Gonstruction of the snapshot space involves solving the local problems for various choices of input 
parameters, and we describe the details below. 
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Coarse mesh 



Figure 2: Illustration of a coarse neighborhood and oversampled domain 


3.1 Snapshot space 

First, we introduce the concept of coarse neighborhood as before. An illustration of coarse neighborhood is 
shown in Figures and respectively, for Darcy and Stokes equations. The snapshot space is composed 
of harmonic extension of fine-grid functions defined on the boundary of Ui excluding the inactive cells in it. 
The local snapshots satisfy boundary conditions imposed on . More precisely, for each fine-grid function, 
which is defined by = Si V/, /c G Jh(oJi), where denotes the fine-grid boundary node on 

For parameter-independent problem, we solve 

inc^AS" (9) 

subject to the boundary condition = S]^{x) on duJi and = 0 on dB^ D cJ^. Note that for the 

du 

differential operator with Neumann boundary condition — = 0, we will use Neumann boundary condition 

on 

instead for the local problems proposed above. 

We denote the space composed of as 

14 nap = span{'*/)®"'‘P : I <l< Mi), 

for each coarse neighborhood where Mi denotes the number of snapshots in the region uji. We emphasize 
that an oversampling strategy (refer to Figure]^ for an illustration of oversampling domain) can be applied 
for the construction of snapshots to obtain a higher accuracy m)- 

Remark 3.1. The snapshots described above are referred to as the harmonic basis in the numerical part. 

Remark 3.2. One can also use all the fine grid functions as snapshots. These snapshots are referred to as 
the spectral basis in the numerical part. 

3.2 Offline space 

This section is devoted to the construction of the offline space via a spectral decompostion. In order to 
construct the offline space Kff, we perform a dimension reduction in the snapshot space using an auxiliary 
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Figure 3: Illustration of a coarse neighborhood for Stokes equations 


spectral decomposition. The main objective is to seek a subspace of the snapshot space such that it can 
approximate any element of the snapshot space in the appropriate sense defined via auxiliary bilinear forms. 
We will consider the following eigenvalue problems in the space of snapshots: 


The definitions of and S°^ for Laplace, elasticity and Stokes equations are listed below. To generate 
the offline space we then choose the smallest Moff eigenvalues from Eqn. ([1Q| ) and form the corresponding 
eigenvectors in the respective snapshot space by setting /c = 1,..., Moff, where 

are the coordinates of the vector We then create the offline matrices 


• For Laplace operator Q, 


Roa = [^f, • • •, V’Moff] • 





and 


^off 



/ snap / snap 
‘r m ‘r n 


( 11 ) 


For elasticity operator (|^, 


= \a°± 


[ (2/x£(V^^^p) : £(C"^P) + AV ■ V ■ V’r^p) , and 

= [c®n] = [ (A+ 

j UJn \ 


( 12 ) 


7 


















o 


o 


o 


o 


o 

o 

o 



o 


o 

o 


Figure 4: Two heterogeneous perforated medium used in the simulations. 


• For Stokes operator 

= [«fn] = / and 

s°^ = [s°J^n]= [ 

Note that for the Stokes operator, multiscale spaces are constructed for velocity field only, and we use 
piecewise constant approximation on the coarse grid for the pressure. 

We remark that some adaptive procedures for choosing the offline basis are proposed and analyzed in 
[IQIIH]. In our future work, we plan to present an analysis of the method which can guide us in choosing the 
local spectral problems. 

4 Simulation results 

In this section, we present simulation results using the framework presented in Sectionj^for Laplace equation, 
elasticity equation and Stokes equations, respectively. We set Q = [0,1] x [0,1] and use two types of perforated 
domains as illustrated in Figure where the perforated regions are circular. Note that we can also 
use perforated regions of other shapes instead and obtained similar results. The computational domain is 
discretized coarsely using uniform triangulation as shown in the left of Figures and where the coarse 
mesh size H = ^. Furthermore, nonuniform triangulation is used inside each coarse triangle element to 
obtain a finer discretization. Examples of this triangulation are displayed on the right of Figures and 
where 2188 and 2364 nonuniform triangle elements are generated, respectively. 

4.1 Laplace equation in perforated domain 

First, we consider the Laplace operator @ imposed with zero Dirichlet boundary condition on the holes 
n dB^ and = 1 on and / = 0. 

The simulation results in perforated domains as shown in Figure are illustrated in Figures and 
respectively. The multiscale solution is obtained in an offline space of dimension 432 (using 12 basis 
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Figure 5: Coarse-scale and fine-scale discretization corresponds to the heterogeneous medium on the left of 
Figure 



Figure 6: Coarse-scale and fine-scale discretization corresponds to the heterogeneous medium on the right 
of Figure 


functions per coarse neighborhood) and the fine-scale reference solution is obtained in a space of dimension 
1187. Compared the fine-scale solution on the left with the coarse-scale solution on the right of the figures, 
we observe that the GMsFEM can approximate the fine-scale solution accurately. 

Besides, two types of snapshot spaces are applied in our simulations with the results shown in Tables 
and|^ The first column of each table shows the number of basis in each coarse node (Nc). The dimensions 
of the offline spaces are given in the second column. The next two columns display the L 2 and Hi relative 
errors for using spectral basis. We can see that the L 2 and Hi relative errors are 0.2% and 2% respectively 
when the dimension of the coarse space is 432. The results applying harmonic extenstion type of basis (Eqn. 

are listed in the last two columns. As displayed in Table the L 2 and Hi relative errors are 1% and 
17% respectively when the dimension of the coarse space is 432. Note that the GMsEEM gives a very good 
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Figure 7: Fine-scale and coarse-scale solution to the Laplace equation in the heterogeneous medium on the 
left of Figure]^ The dimension of the coarse space is 432. 
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Figure 8: Fine-scale and coarse-scale solution to the Laplace equation in the heterogeneous medium on the 
right of Figure The dimension of the coarse space is 432. 


solution in both cases with only about 36% of unknown compared to the fine-scale solver. We also observe 
in Tables [2 and that a rapid decay in the errors as we increase the number of basis functions. 

4.2 Elasticity equation in perforated domain 

Next, we consider elasticity operator <§• We use zero displacements = 0 on the holes, = 0, = 0 
on the left boundary, da, = 0, = 0 on the bottom boundary and da, = 0, = 0 on the top and right 
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Nc 

dim 

harmonic basis 
L 2 Hi 

spectral basis 
L2 Hi 

1 

36 

0.18 

0.75 

0.19 

0.63 

2 

72 

0.12 

0.65 

0.10 

0.45 

4 

144 

0.06 

0.48 

0.02 

0.18 

8 

288 

0.03 

0.35 

0.01 

0.09 

12 

432 

0.01 

0.17 

0.002 

0.02 

16 

576 

0.004 

0.04 

0.001 

0.009 


Table 1: Numerical tests for Laplace operator in heterogeneous medium shown on the left of Figure 
Fine-scale problem dimension is 1187. 


Nc 

dim 

harmonic basis 
L2 Hi 

spectral basis 
L2 Hi 

1 

36 

0.10 

0.89 

0.06 

0.34 

2 

72 

0.09 

0.88 

0.03 

0.25 

4 

144 

0.05 

0.65 

0.01 

0.10 

8 

288 

0.03 

0.52 

0.005 

0.05 

12 

432 

0.01 

0.29 

0.001 

0.02 

16 

576 

0.003 

0.06 

0.0006 

0.009 


Table 2: Numerical tests for Laplace operator in heterogeneous medium shown on the right of Figure 
Fine-scale problem dimension is 1269. 


boundaries. Here, u = {ux^Uy) and a = {(Tx^cfy)^ The source term is defined by / = (10^, 10^), the elastic 
modulus is given hy E = 10^, Poisson’s ratio is = 0.22, where 

E Ev 

^ 2(1+ (1 + z/)(l - 2z/) ‘ 

The fine-scale solution and coarse-scale solution corresponding to different perforated domains in Figure 
are presented in Figures [9| and pT| The fine-scale displacement is displayed on the left, while the coarse-scale 
displacement is displayed on the right. The multiscale solution is obtained in an offline space of dimension 
864 (using 12 basis functions per coarse neighborhood for each component of the displacement) and the 
fine-scale reference solution is obtained in a space of dimension 2374. Comparing the fine-scale solution with 
the coarse-scale solution in Figures and we can observe a good accuracy. 

Furthermore, the relative errors for different dimensions of coarse spaces with harmonic basis functions 
(Eqn. and spectral basis functions are shown in Tables andrespectively. As displayed in Table 
the 1/2 and Hi relative errors using harmonic basis are 3% and 14%, respectively, when the dimension of 
the coarse space is 1152 for the perforated domain shown in the left of Figure In addition, for the same 
dimension of the coarse space and the same domain, the L 2 and Hi relative errors using spectral basis are 
0.5% and 5%, respectively Note that the dimension of the fine-scale solution is 2374. Thus, we see that the 
proposed GMsFEM gives accurate solutions. We also observe that GMsEEM gives more accurate solution 
when more basis are included in the offline space. We remark that the results (see Table follow a similar 
pattern for the domain shown on the right of Eigurej^ 


4.3 Stokes equations in perforated domain 

In our final example, we consider the Stokes operator with zero velocity u = (0, 0) on D dB^ and 
u = (1,0) on 90. The numerical results corresponding to different perforated domains in Eigurej^are listed 
in EigurespT|and 12 The fine-scale solution and coarse-scale solution are depicted on the left and right of 
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Figure 9: The fine-scale and coarse-scale solutions of the elasticity equation correspond to the heterogeneous 
perforated domain shown on the left of Figure]^ The dimension of the fine-scale solution is 2374 and the 
dimension of the coarse space is 864. 


Ne 

dim 

harmonic basis 
L2 Hi 

spectral basis 
L2 H, 

1 

72 

0.65 

0.74 

0.58 

0.75 

2 

144 

0.41 

0.57 

0.21 

0.46 

4 

288 

0.32 

0.48 

0.07 

0.26 

8 

576 

0.16 

0.32 

0.02 

0.13 

12 

864 

0.11 

0.25 

0.009 

0.08 

16 

1152 

0.03 

0.14 

0.005 

0.05 

20 

1440 

0.01 

0.04 

0.003 

0.03 


Table 3: Numerical tests for elasticity equations in heterogeneous medium shown on the left of Figure 
The fine-scale problem dimension is 2374. 
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Figure 10: The fine-scale and coarse-scale solutions of the elasticity equation correspond to the heterogeneous 
perforated domain shown on the right of Figure The dimension of the fine-scale solution is 2538 and the 
dimension of the coarse space is 900. 


Nc 

dim 

harmonic basis 
L2 Hi 

spectral basis 
L2 Hi 

1 

72 

0.79 

0.86 

0.65 

0.81 

2 

144 

0.53 

0.69 

0.11 

0.35 

4 

288 

0.45 

0.62 

0.04 

0.20 

8 

576 

0.22 

0.41 

0.01 

0.12 

12 

864 

0.16 

0.34 

0.007 

0.07 

16 

1152 

0.06 

0.20 

0.003 

0.05 

20 

1440 

0.01 

0.12 

0.002 

0.04 


Table 4: Numerical tests for elasticity equations in heterogeneous medium shown on the right of Figure 
The fine-scale problem dimension is 2538. 
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Figures pT] and \T^ respectively. To improve the accuracy, we have enriched pressure spaces by considering 
a splitting algorithm [5]. Our preliminary numerical results show an improvement. 


Nc 

dim 

harmonic basis 
L 2 H, 

I 

108 

0.85 

0.88 

2 

180 

0.68 

0.78 

4 

324 

0.40 

0.69 

8 

612 

0.26 

0.67 

12 

900 

0.23 

0.60 

16 

1188 

0.23 

0.58 


Table 5: Numerical tests for Stokes operator in heterogeneous medium shown on the left of Figure Fine- 
scale problem dimension is 3561. 


Nc 

dim 

harmonic basis 
L 2 H, 

1 

108 

0.65 

0.87 

2 

180 

0.56 

0.82 

4 

324 

0.31 

0.71 

8 

612 

0.16 

0.62 

12 

900 

0.13 

0.46 

16 

1188 

0.13 

0.44 


Table 6: Numerical tests for Stokes operator in heterogeneous medium shown on the right of Figure 
Fine-scale problem dimension is 3807. 

The 1/2 and Hi relative errors in different perforated domains in Figure are listed in Tables and 
respectively. We observe that the L 2 and Hi relative errors are 23% and 58%, respectively. As more basis 
are taken in the construction of the offline space, the L 2 and Hi relative errors become smaller. In general, 
enriching pressure spaces and constructing associated velocity spaces can improve the accuracy. This is 
currently under investigation. 

4.4 Randomized snapshots for GMsFEM 

In this subsection, we will investigate the oversampling randomized algorithm proposed in [2] (shown in Table 
[^. The advantage of this algorithm lies in the fact that a much fewer number of snapshot basis functions are 
calculated maintaining a good accuracy to the solution space concurrently. Besides, the oversampling strategy 
is used to reduce the mismatching effects of boundary conditions imposed artificially in the construction of 
snapshot basis functions. For the sake of brevity, the results for elasticity equation are shown only. 

The simulation results are presented in Tables and for elasticity problem. The results without 
oversampling are shown on the top of each table, while those with the oversampling strategy are shown on 
the bottom part of each table. In our simulation, we set the oversampling size t = 2 (i.e., two extra fine-grid 
blocks are added to the original coarse region) and a buffer number = 4 for each coarse neighborhood 
uji^ i.e., we will generate 4 more snapshot basis functions in each coarse neighborhood. In the numerical 
results, we report the fraction of the snapshots computed compared to all snapshot vectors. For the sake of 
completeness, we list the algorithm in Table 

Comparing the results on the left columns (using full snapshot space) of Tables and with those on 
the right columns that correspond to the randomized snapshots, we observe that the randomized algorithm 
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Figure 11: The fine-scale and coarse-scale solutions of the pressure, x-component and y-component of the 
velocity correspond to the heterogeneous perforated domain shown on the left of Figure The dimension 
of the fine-scale solution is 3561. 
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Figure 12: The fine-scale and coarse-scale solutions of the pressure, x-component and y-component of the 
velocity correspond to the heterogeneous perforated domain shown on the right of Figure The dimension 
of the fine-scale solution is 3807. 
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Table 7: Randomized GMsFEM Algorithm 


Input: 

output: 

1 . 

2 . 

3. 

4. 


Fine grid size h, coarse grid size i7, oversampling size t, buffer number for each 
the number of local basis functions for each 
Coarse-scale solution uh- 

Generate oversampling region for each coarse block: T^, T^, and uj ^; 

Generate random vectors ri and obtain randomized snapshots in ; 

Add a snapshot that represents the constant function on ; 

Obtain k^l offline basis by a spectral decomposition restricted to the original region; 
Gonstruct multiscale basis functions and solve it. 


Ne 

full snapshots 
L2 Hi 

randomized snapshots 
L2 Hi 

without oversampling, Wi 


100 % 

39.7 % 

1 

0.572 

0.753 

0.764 

0.853 

2 

0.217 

0.466 

0.592 

0.742 

4 

0.071 

0.261 

0.317 

0.529 

8 

0.023 

0.136 

0.173 

0.386 

12 

0.009 

0.079 

0.101 

0.286 

16 

0.005 

0.054 

0.055 

0.214 


with oversampling, 

= Wi + 2 


100 % 

25.2 % 

1 

0.506 

0.716 

0.533 

0.729 

2 

0.223 

0.472 

0.203 

0.450 

4 

0.069 

0.258 

0.065 

0.252 

8 

0.025 

0.143 

0.025 

0.146 

12 

0.010 

0.083 

0.013 

0.095 

16 

0.006 

0.058 

0.007 

0.069 


Table 8: Numerical tests for elasticity in heterogeneous media shown on the left of Figure Randomized 
oversampling for GMsFEM. 


provides a nearly similar results, particularly with oversampling strategy. This is consistent with our previous 
observations for heterogeneous problems. Oversampling avoids the oscillations on the boundary conditions 
due to the randomness. The effect of oversampling strategy is much more obvious for the randomized 
snapshot space. 

5 Conclusion 

In this paper, we develop a generalized multiscale finite element framework for problems in perforated 
domains. Our approach follows GMsFEM. The main contributions of this paper are the development of 
snapshot space and local spectral problems for solving problems in perforated domains with multiple scales 
and no scale separation. Our approaches differ from previously developed GMsFEM techniques for hetero¬ 
geneous problems. In particular, the snapshot vectors and local eigenvalue problems need to be developed 
in a multiscale domain and take into account the boundaries that may be disconnected. We show that using 
GMsFEM framework, we can propose a unified framework for solving problems in perforated domains. In 
the paper, we discuss three applications: (1) Laplace equation in perforated domain; (2) elasticity equation 
in perforated domain; and (3) Stokes equations in perforated domain. We present some preliminary numer- 
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N 

full snapshots 

randomized snapshots 


L2 

Hi 

L2 

Hi 


without oversampling 



100 % 

39.2 % 

1 

0.692 

0.822 

0.858 

0.908 

2 

0.116 

0.351 

0.673 

0.798 

4 

0.039 

0.207 

0.510 

0.678 

8 

0.016 

0.128 

0.358 

0.551 

12 

0.007 

0.078 

0.204 

0.418 

16 

0.003 

0.055 

0.132 

0.334 


with oversampling, = 

Wi ^2 


100 % 

23.3 % 

1 

0.665 

0.800 

0.611 

0.762 

2 

0.115 

0.355 

0.131 

0.383 

4 

0.038 

0.214 

0.049 

0.241 

8 

0.015 

0.135 

0.020 

0.157 

12 

0.006 

0.087 

0.008 

0.108 

16 

0.006 

0.089 

0.008 

0.102 


Table 9: Numerical tests for elasticity in heterogeneous media shown on the right of Figure Randomized 
oversampling for GMsFEM. 

ical results that show that one can efficiently solve these problems on a coarse grid using fewer degrees of 
freedom. We also discuss the use of randomized snapshots to reduce the offline computational cost associated 
with computing the snapshot space. In our future work, we plan to present analysis and design new more 
efficient coarse spaces based on the analysis. 
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